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We study the behavior of a Bose-Einstein condensate held in an optical lattice. We first show how a self- 
trapping transition can be induced in the system by either increasing the number of atoms occupying a lattice 
site, or by raising the interaction strength above a critical value. We then investigate how applying a periodic 
driving potential to the self-trapped state can be used to coherently control the emission of a precise number of 
correlated bosons from the trapping- site. This allows the formation and transport of entangled bosonic states, 
which are of great relevance to novel technologies such as quantum information processing. 
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Introduction Ultracold atoms held in optical lattices are cur- 
rently at the center of intense theoretical and experimental in- 
vestigation. The experimental parameters of these systems 
can be controlled extremely precisely, and in addition, the 
high degree of isolation from the environment permits their 
quantum dynamics to remain coherent over long timescales. 
Consequently such systems provide an attractive way of in- 
vestigating quantum many-body physics, and it has been pro- 
posed |1] to use them to simulate strongly correlated quan- 
tum systems of interest in other areas of physics such as high- 
temperature superconductivity |2], They also provide an ex- 
cellent starting point for engineering and manipulating entan- 
gled states, which are vital for implementations of quantum 
information processing. 

Bosons confined in an optical lattice provide an almost ideal 
realization |3] of the Bose-Hubbard (BH) model, and a recent 
pioneering experiment directly observed the quantum phase 
transition between a superfluid and a Mott insulator |4]. The- 
oretical work |5j,^ has shown how this transition may be in- 
duced in an alternative way: by applying an additional os- 
cillatory driving field to suppress the inter-site tunneling by 
means of the quantum interference effect termed coherent de- 
struction of tunneling (CDT) |7]. In this work we show how 
such an oscillating driving field can be used to manipulate the 
dynamics of a different ground-state of the BH model - the 
self-trapped state. In particular we demonstrate an effect anal- 
ogous to photon-assisted tunneling, in which certain driving 
frequencies induce a coherent oscillation of an integer number 
of bosons between the trapping-site and its nearest-neighbors. 
Combining this effect with CDT allows us to control the emis- 
sion of a definite number of bosons from the trapping site and 
manipulate their speed of propagation through the optical lat- 
tice, and thus enables the self-trapped state to be used as a 
quantum beam-splitter or as a coherent source of entangled 
bosons. 

Static properties The BH model is described by the Hamil- 
tonian 



(jM) 



where ajja'j are the standard annihilation/creation operators 
for a boson on site j, nj - a\aj is the number operator, J is 
the tunneling amplitude between neighboring sites, and U is 



the repulsion between a pair of bosons occupying the same 
site. Here we initialize the system in a state in which all the 
bosons occupy a single lattice site. As the Hubbard interac- 
tion is repulsive, it might be thought that such a state would 
be extremely unstable. Surprisingly, however, this is not nec- 
essarily the case and, depending on the strength of the inter- 
action and the filling of the lattice site, this highly-localized 
configuration is able to persist for long times. In such cases 
the bosons are said to be "self-trapped" 0]. 

Self-trapping has been recently observed experimentally in 
Bose-Einstein condensates of roughly 1000 atoms LS^ J^O], and 
can be understood quaUtatively by an energetics argument. 
The presence of the optical lattice causes the energy spectrum 
of non-interacting bosons to be confined to a Bloch band of 
width 47. Consequently, if the potential energy per particle 
of the trapped condensate is much higher than this, it cannot 
be converted into the kinetic energy of free bosons and the 
trapped state cannot decay. Its stability thus depends critically 
on the absence of dissipative processes in optical lattice sys- 
tems. 

To describe this effect quantitatively, we consider a 2-site 
BH model holding bosons. If the system is initialized in 
the state lA', 0) (that is, a Fock state with bosons occupying 
one site with the other site empty), then the primary tunneling 
process will be to the state |A^- 1, 1). If we truncate the Hilbert 
space to just these two states, we obtain an effective two-level 
model 



ViN) J y/N 
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where V(n) - Un{n - l)/2 is the potential energy of n bosons 
occupying one site. It is useful to visualize the time evolu- 
tion of the system geometrically by making use of the Bloch 
sphere representation. Parameterizing 7^2-lev ^'^ terms of the 
Pauli matrices 

U(N - If r- U(N - 1) 

^2-lev = 2 ^ + 2 '^'^ 

reveals that we can interpret it as an interaction between 
the Bloch vector and a fictitious magnetic field B = 
{J y/N, 0, U(N- l)/2). Thus under the influence of the Hamil- 
tonian the Bloch vector will simply make a Larmor orbit about 
B. This form of time evolution is shown in Fig^ for a 



2 



(a) 

^'^ 





FIG. 1: Time evolution of a two-site system containing 7 bosons. 
The system is initialized in the state |7, 0) which corresponds to the 
north pole of the Bloch sphere, and evolves under the full BH Hamil- 
tonian 0. To the left we show the evolution on the Bloch sphere, 
projected onto the (.s\, Sy) plane, and to the right the overlap of the 
system with the initial state (a) Strong interaction, U/J = 8. 

The system is described well by the two-level model J3j, and peri- 
odically traces out a circle on the Bloch sphere. Correspondingly, p, 
makes high-frequency oscillations of small amplitude, and never at- 
tains a value smaller than 0.99, indicating that self-trapping is occur- 
ring, (b) Intermediate interaction U/J = 4. The radius of the circular 
orbit is larger, but the two-level approximation is still good. The os- 
cillations in pj are now of larger amplitude showing that the trapping 
effect is reduced, (c) Weak interaction, U/J = 1. The two-level 
approximation now breaks down, and the system's time-evolution is 
correspondingly more complicated. The value of rapidly drops as 
the trapping effect has now been completely lost. 



strongly interacting (U - 87) 7-boson system. It can be 
clearly seen that, as expected, the Bloch vector traces out a 
periodic circular orbit centered on B. Since for these parame- 
ters B is almost parallel to the s. axis, the radius of this orbit 
is rather small, and thus the system exhibits a high degree of 
self-trapping. 

To assess the degree of self-trapping more precisely we 
measure the overlap of the system's state with the initial state, 
Piit) - K^AiliACO)!^. since when self-trapping occurs this quan- 
tity is restricted to values close to unity. The circular motion 
shown in Fig[^ equates to a very small-amplitude sinusoidal 
oscillation of /?,. As U is reduced, the angle between B and the 
s, axis increases, and as a result the radius of the Larmor orbit 
made by the Bloch vector increases (FigOJ)). Consequently 
the degree of trapping is reduced, and the amplitude of oscilla- 
tions of Pi is larger If U is decreased even further (Fig[lJ;) the 
self-trapping effect is lost, and the two-level approximation 
breaks down. In this case the Bloch vector rapidly leaves the 
surface of the Bloch sphere, and its erratic time-evolution cor- 
responds to an irregular quasi-periodic behavior of which 
can take very low values. 



In the inset of Fig|2we plot the minimum value of at- 
tained in the 7-boson system as a function of the interaction 
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FIG. 2: The dashed-line plots the value of the Hubbard-interaction, 
Uc, at which the dynamics of an undriven two-site system crosses 
over from the untrapped to the self-trapped regime (Eq|3J. For large 
particle numbers Uc ~ 1/ V^, but rises sharply as the number of 
bosons is reduced. Inset: Comparison of the two-level approximation 
with the full BH Hamiltonian for a 7-boson system. For U/J< I 
the 2-level approximation breaks down, but becomes increasingly 
accurate as U is increased. 



strength. For small values of U no self-trapping occurs, and 
thus takes a value of zero. As U is increased, however, 
the oscillations in pi are quenched, and for U/J > 3 it can be 
seen that the effective two-level model describes the dynamics 
extremely well. From Eql^we can obtain a criterion for the 
crossover to the self-trapped regime, by defining the transition 
to occur when p™" drops below a value a. This yields a value 
for the critical value of U, 
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J 
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This boundary is plotted in Fig|2for a = 0.99 - that is, when 
less than 1% of the initial state is able to leak to the other site. 
Although the self-trapping regime is more easily achievable 
for large boson numbers since Uc ~ 1/ ^/N, it is in princi- 
ple possible in a system of just two bosons, if the interac- 
tion strength can be raised sufficiently high. This has recently 
been achieved experimentally 111 IB in a gas of trapped rubid- 
ium atoms. 



Dynamical properties We now consider dynamically con- 
trolling the self-trapped state by applying a harmonic driving 
potential 



H(t) - Hbh + K sin cot ^ j nj. 



(5) 



where K is the amplitude of the driving field, and w is its fre- 
quency. We first consider the case of extracting a single bo- 
son from the trapping site. As the occupation of the trap site 
changes from N to N-l there is a corresponding loss of poten- 
tial energy AE - Vn - VOv-i = U(N - 1), where V„ is defined 
as in Eq|2] When a system possesses such a large energy gap. 
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FIG. 3: Response of a 5-site system, holding 5 bosons (U = 16), to 
a sinusoidal driving field, (a) For a driving frequency to = 4(7, the 
value of p™"^ rapidly drops as K is increased from zero, indicating 
that the initial state is rapidly destroyed by photon-assisted tunneling. 
At a sequence of sharp peaks centered on the zeros of (K/a)) (ver- 
tical dashed lines) coherent destruction of tunneling (CDT) instead 
causes the system to be frozen in its initial state, (b) Time depen- 
dence of the occupation numbers of the lattice sites. For t < 10, 
K/(i} is set to be the first zero of , and CDT suppresses any oscil- 
lations. For t > 10 we set K/<i> = 2.40 which produces a sinusoidal 
oscillation of a single boson from the central site (n = 0) to its two 
neighbors (n = ±1). (c) As in (a) but for a frequency of io = 3U. 
Again p™"^ peaks at the zeros of Ji^K/o)), although the peaks are 
rather broader, (d) As in (b), but for to = 3U. When K/to switches 
to a value of 2.40 the number densities show a slow sinusoidal oscil- 
lation, in which the occupation of the central site varies between five 
and three, i.e. this driving frequency induces a two-boson oscillation. 



Floquet analysis may be used to show p6|,|12l] that extremely 
fine control over the tunneling dynamics is possible at multi- 
photon resonances, that is, when ma> — AE, m — 1,2 . . .. In 
general when this condition is satisfied, the system is able to 
exchange energy with the driving field to overcome the energy 
gap, and so tunneling is restored 1 13]. However, at particular 
values of the amplitude of the driving field, CDT will occur 
when the Floquet quasienergies of the system become degen- 
erate, and correspondingly the dynamics of the system will 
be frozen. For a sinusoidal driving potential these degenera- 
cies occur at the zeros of J',n(K/(D), the mth Bessel function 
of the first kind. Thus at a photon resonance it is possible to 
produce dramatic differences in the tunneling rate by making 
small changes in the amplitude of the driving field to move the 
system between CDT and photon assisted tunneling. 

In Figl^ we plot obtained in a 5-boson system 

driven at a frequency of w = 4U. This corresponds to the 
m - 1 photon resonance. For K - Q the system is self- 
trapped, and consequently p, remains near unity. Increas- 
ing K, however, causes the value of p™" to rapidly drop 
to zero, demonstrating how photon assisted tunneling over- 
comes the self-trapping effect. As K is increased further, 
^mm exjjj^ifs ^ number of extremely sharp peaks centered on 
K/a> = 3.83, 7.01 and 10.17 - the zeros ofJi(K/aj). 



Away from these zeros, the driving field causes a single 
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FIG. 4: (a) Schematic version of the process illustrated in Figlsjj. 
The driving field excites a boson from the trapping site (n = 0). If 
Jo(K/tL)) = 0, CDT suppresses the single boson tunneling processes 
(marked by X), and so a Rabi-like oscillation occurs between the 
trapping site and its neighbors, (b) Changing the driving parame- 
ters to stabilize the self-trapping isolates the remaining bosons in the 
trapping site, but the escaped particle is able to propagate through the 
optical lattice with a renomalized tunneling amplitude J^f^. 



boson to diffuse symmetrically from the trapping site to its 
neighboring lattice sites, and from there to continue propa- 
gating through the optical lattice with a renormalized tunnel- 
ing 1 7] Jgfi - JJ'oiKIco). If we therefore choose a value of 
K such that J^^ - 0, the particle will not be able to prop- 
agate further, and will thus just make a Rabi-like oscillation 
between the trapping site and its neighbors. This situation 
is illustrated schematically in Fig|^. We show in Fig|3j5 the 
time dependence of the occupation of the trapping site and its 
neighbor Initially we set K/cj - 3.83 (a zero of J'l), and no 
tunneling occurs: for this value of K/a> CDT reinforces the 
self-trapping effect. At f = 10 we alter the amplitude of the 
driving to K/co = 2.40 (a zero of J'q) and a very clear parti- 
cle oscillation occurs, in which the occupation of the trapping 
site cycles between 5 and 4, and that of the neighboring sites 
moves between and 0.5. In this sense the optical lattice is 
acting as a particle beam-splitter, dividing a single boson into 
a superposition of left and right-moving components. 



We have so far considered the case of extracting a single 
boson from the trapping site. We can however apply a similar 
method to induce the emission of an integer number of bosons, 
Nem, by noting that the energy difference per particle takes the 
remarkably simple form 



AE/N,„ = [VN-{VM-N„, + VNj]/Ne. 
= U{N-Ne,n) 



(6) 



The self-trapped state thus contains a ladder of equally-spaced 
energy levels, separated by U, and so by driving the system at 
the correct frequency we can induce the emission of a given 
number of particles. In Fig|3l; we show the response of the 5- 
boson system to a driving field of frequency u) - 3U, thereby 
inducing the emission of two bosons. As before we can ob- 
serve peaks in p™" centered on the zeros of J'i{K/cj) at 
which CDT occurs, while between them photon-assisted tun- 
neling causes p™"^ to take low values. FiglSjl shows the effect 
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FIG. 5: Above: Time evolution of the particle density of the trapping 
site (n = 0) and its neighbor {n = 1) for a 5-boson system. The 
driving parameters are chosen (see text) to produce the emission of 
single bosons at f 5 and 25, and thus the occupation of the trapping 
site falls in steps as 5 — > 4 — > 3. Below: Particle density of the 
system. The emitted particle bursts move away from the central site 
at a roughly constant speed, and disperse slightly as they propagate 
through the lattice. 

of switching the amplitude of the driving field to a value of 
Kju) - 2.40. We can again see that this has the effect of in- 
ducing a Rabi-like oscillation between the trapping site and 
its neighbors, but in this case the oscillation indeed consists 
of two bosons, and has a longer period. 

In Fig|5]we show how combining the PAT effect with CDT 
allows the population of the self-trapped state to be reduced 
step-by-step. The system is initialized with 5 bosons in the 



central site, and is initially driven at w = AU. Driving the sys- 
tem at Kjto - 3.80 induces CDT which suppresses the num- 
ber variance of the trapping site and stabilizes the self-trapped 
state. At f = 5 the value of Kju is changed to 2.40 which in- 
duces the Rabi-type oscillation between the trapping site and 
its neighbors, as schematically shown in Fig|^. After a half- 
integer number of these oscillations we then alter the driving 
parameters to w = 3U, Kjio - 3.80 which traps the remain- 
ing 4 bosons in the trapping site. For these parameters, how- 
ever, the single-particle tunneling is not quenched and so the 
ejected boson is able to propagate through the optical lattice 
away from the trapping site, as shown in Fig|^. The emitted 
particle smears out to an extent as it moves through the lattice 
(see Fig|5] lower panel), but nonetheless the atom-pulse re- 
mains quite clearly defined after propagating through several 
lattice spacings. By repeating this procedure with appropriate 
driving frquencies we can successively reduce the occupation 
of the trapping site in integer steps, and thereby produce a se- 
quence of well-defined, phase-coherent atom-pulses. 



Conclusions We have shown how self-trapping arises in a 
Bose-Einstein condensate confined to an optical lattice. Ap- 
plying a resonant driving field to the self-trapped state can ei- 
ther stabilize the trapping (when CDT occurs), or can induce a 
Rabi-like particle oscillation. The interplay between these ef- 
fects makes it possible to control the emission and propagation 
of a precisely-defined number of particles, and thus enables 
the self-trapped state to be used as a particle beam-splitter or 
a source of mesoscopic entangled states, which have many 
possible appUcations in quantum information. 
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